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' Abstract 

We summarize recent results obtained in the Dyson-Schwinger formalism to study the chiral 
and deconfinement phase transitions of quenched and unquenched QCD at finite temperature and 
^ I chemical potential. In the quenched case we compare SU{2) and SU{3) gauge theories by taking 

lattice data for the gluon as an input for the quark Dyson-Schwinger equation. As compared to 
previous investigations we find a clearer distinction between the second order transition of the 
. two-color theory and the (weak) first order transition of the three-color gauge theory. We then 

QQ I extend this study to unquenched QCD at finite chemical potential by taking matter effects to the 

gluon into account and investigate the order of the chiral phase transition and the behavior of the 
deconfinement transition. What we find are coinciding phase transitions up to a critical endpoint 
which is located at large chemical potential. 



1 Introduction 



The behavior of quantum chromodynamics at large temperatures and densities received a lot of attention 
■ over the past years and is an ongoing research program from both, theoretical and experimental side. 
At vanishing chemical potential lattice QCD has shown the existence of a crossover from the chiral 
symmetry broken and confined hadronic phase to the phase of the (approximately) chiral symmetric 
and deconfined quark-gluon plasma. While this has been confirmed many times, the behavior at finite 
chemical potential is still under intense debate. Lattice methods have a limited applicability here due 
to the fermion sign problem. 

Models like the Nambu-Jona-Lasinio (NJL) and the quark-meson (QM) model have so far been the 
main source of studies at large chemical potential, and established a scenario where the chiral crossover 
turns into a first order phase transition at a critical endpoint. In Polyakov loop extended versions these 
models (PNJL and PQM) have also been used to study the confinement /deconfinement transition, 
which may or may not coincide with the chiral transition [U El |3l H]. At large chemical potentials 
and relatively small temperatures there may also exist some new phases, e.g. color superconductors, 
inhomogeneous O [6] or quarkyonic [7J phases. 

Another direct approach to non-perturbative QCD without the sign problem is the framework of 
Dyson-Schwinger equations (DSEs) [H El [1^ and the functional renormalisation group [HI [12]. QCD 
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with two degenerate quark flavors has been studies in Ref. [T0| by solving the coupled system of quark 
and gluon DSEs using quenched lattice data for the gluon propagator as input. Within this truncation 
scheme the behavior of the chiral and deconfinement transitions at finite chemical potential have been 
investigated using the first calculation of the dressed Polyakov loop in this region of the QCD phase 
diagram. In this proceedings contribution we give an overview of the employed truncation scheme and 
summarize the corresponding results. 



2 Order parameters for chiral symmetry breaking and con- 
finement 

The central object of our investigations is the in-medium quark propagator which can be decomposed 
as 

where n is the quark chemical potential and a;„ = 7rT(2n + l) are the Matsubara modes in the imaginary 
time formalism with temperature T. The functions A, C and B dress the vector and scalar part of the 
propagator which we calculate from the corresponding DSE. The bare propagator 5*0 at quark mass 
m is characterized hj A = C = 1 and B = m. A quark propagator with a non-vanishing B function 
corresponds to a phase of broken chiral symmetry, either explicitly by a bare quark mass or dynamically 
by quantum effects. 

A possible order parameter for chiral symmetry breaking is the quark condensate 

which can be calculated directly from a solution of the quark DSE. The condensate in this definition 
is divergent with mA^ and m^A, but since these terms do not depend on temperature and chemical 
potential the condensate can still be used as an order parameter. We work with approximately physical 
quark masses and therefore expect to find a crossover at small chemical potentials. To define the 
pseudo-critical temperature in this case we use the susceptibility 

^(^V") /ON 

X = -j^ , (3) 

am 

and determine its maximum to find Tg. The divergent terms in the condensate only lead to an offset 
in X, without changing its maximum. The quark condensate in the chiral limit has been determined in 
Ref. [H], where critical scaling beyond the mean field level at the second order chiral phase transition 
has been studied. 

Constructing an order parameter for confinement that is accessible with functional methods is a 
more challenging task. In the quenched case the Polyakov loop is an accepted order parameter which 
is known to show a (weak) first order phase transition for QCD at a temperature of approximately 
270 MeV and a second order transition for the two-color gauge theory at approximately 300 MeV. 
The Polyakov-loop expectation value is sensitive to breaking of centre symmetry, which occurs in the 
deconfined phase, while the confined phase is centre symmetric. However, if finite quark masses are 
taken into account within QCD, centre symmetry is broken explicitly (just like chiral symmetry). There 
are, however, interesting indications that this might not be the case if the larger centre symmetry of 
the full standard model is taken care of [13j. 

In [151 [13 [H] the so-called dual condensates have been proposed as order parameters for centre 
symmetry breaking. They are defined as 
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Figure 1: The Dyson-Schwinger equation for the quark propagator. Dots denote dressed objects. 




Figure 2: The truncated gluon DSE. The black dot denotes the full, unquenched propagator, while the 
grey dot denotes the quenched propagator. 

s„= r%-'^-m^, (4) 

Jo 2,71 

where (p G [0, 27r[ is a parameter for [/(l)-valued boundary conditions of the quark fields: ijj{x, 1/T) = 
e^'-'^iplx, 0). The physical boundary condition is = tt. The quantity S„ contains all loops of connections 
winding n times around the Euclidean time direction. For n = ±1 this also contains the Polyakov loop, 
and E-ti has therefore been named the 'dressed Polyakov loop'. It also contains all kinds of loops which 
are not straight in the time direction but contain detours, but these loops are 1/m suppressed. In the 
m — )■ oo limit the normal Polyakov loop is recovered. At finite chemical potential S+i and S_i are 
not equal and correspond to the dressed Polyakov loop and its conjugate. Since the dressed Polyakov 
loop is sensitive to centre symmetry breaking, it allows us to calculate the confinement /deconfinement 
transition from solutions of the quark DSE without having to deal with an ansatz for a Polyakov-loop 
potential as necessary in the PNJL and PQM models. 

Since the deconfinement transition will be a crossover in the unquenched case, we use the maximum 
of to define the pseudo-critical temperature, similar to the chiral transition. 



3 Dyson-Schwinger equations 

Fig. [1] displays the DSE for the quark propagator. The quark self-energy depends on the fully dressed 
gluon and quark-gluon vertex, which we need to specify in order to solve the equation self-consistently. 
For the in-medium gluon propagator we take two steps, first we investigate quenched QCD where lattice 
calculations are up to now the most reliable source for the temperature dependent gluon propagator, 
and then we will introduce unquenching effects by resorting to the gluon DSE. 

The unquenched gluon DSE in the truncation we use here is depicted in Fig. O The full propagator 
is given by the lattice results for the quenched propagator and the quark loop. This is only an ap- 
proximation, since all diagrams where quark loops appear inside the Yang-Mills part of the self-energy 
are neglected. In vacuum this leads to an over-estimation of the quark loop on the few percent level. 
Assuming this still holds at finite temperature, we expect to underestimate the critical temperatures 
by up to 10 MeV. 

The quark loop shown in Fig. [2] is given by 

^AV) = E / 7^ Tr [S{q)g^,S{k)gV^] , (5) 
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where q = p+k, Z\p is the vertex renormahzation factor and Y^, the full quark-gluon vertex. In principle 
the system of quark and gluon DSEs can now be solved self-consistently but as a first approximation we 
will treat the quark loop semi-perturbatively by taking bare quarks but a dressed vertex. This allows 
us to use the hard-thermal loop expression multiplied by the vertex dressing function, defined below. 
The HTL approximation is well justified above the critical temperature where quark dressing effects are 
small, but needs to be corrected in the chiral broken phase. A calculation with a fully dressed quark 
loop is work in progress. 

Finally we have to specify our choice of the quark-gluon vertex, which appears in the quark self- 
energy and in the quark loop. We use the same construction as in [9]. It is given by 



C{v) + C{q) 



+ 5, 



A{v)\A{qy 



11,1 



di 



'/3oa{i2) ln[gVA2 + l]' 
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An 



(6) 
(7) 



In the UV resummed perturbation theory is recovered with f3o = (llA'c — 2A'/)/3, = 0.3, 6 = 
—9Nc/ (44iVc — 8A^/) and the scale A = 1.4 GeV, which is inherited from the lattice gluon. In the IR we 
have two parameters, di and d2, which have been fixed to d2 = 0.5 GeV^ and di = 7.6 GeV^ for SU{2) 
and di = 4.6 GeV^ for SU{3). 



4 Results 



4.1 Quenched QCD 

As already mentioned above, the phase transition of quenched QCD (i.e. without the quark-loop 
contributions in the Yang-Mills sector) can be determined from the quark DSE using quenched lattice 
data for the temperature dependent gluon propagator as input. Of course, the quality of the results 
then depends on the statistic and systematic error of the lattice. Starting from the pioneering work of 
Ref. [18], these have been improved in and analyzed in more detail in recent works [191 [20l [211 [22] . 
general, however, it seems fair to say that in particular systematic errors due to volume and discretization 
artifacts at small momenta are not yet well under control. This is particularly true in the vicinity of the 
critical temperature. Consequently, it proved difficult to distinguish the order of the phase transition 
between the two-color and three-color cases investigated in Ref. [5]. Here we present updated and 
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Figure 3: Dressed Polyakov loop and quark condensate for SU{2) (left) and SU{3) (right). 
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Figure 4: Dressed Polyakov loop and quark condensate in two flavor QCD as function of temperature 
at zero chemical potential. 

improved results for the quark condensate and the dressed Polyakov loop using the high-statistics 
lattice data of Ref. [22] as input, which are also carried out on a much finer temperature grid than the 
ones of [9J. 

Fig. |3] shows how the dressed Polyakov loop and the quark condensate change with temperature in 
the cases of two and three colors. For the normalization of the temperature scale we use the transition 
temperatures which have been determined from the Polyakov loop on the lattice and compare with 
our results obtained from the quark DSE. Indeed, both order parameters show a rapid change at Tc, 
signaling the (approximate) restoration of chiral symmetry and breaking of centre symmetry at the very 
same temperature in agreement with the critical temperature determined from the lattice. With the 
finer temperature grid and the better statistics compared to [9] , the behavior of the dressed Polyakov 
loop is now clearly distinguishable between the SU{2) and SU{3) cases, pointing towards a second order 
phase transition for SU{2) and a weak first order for SU{3), again in agreement with the expectations. 
The situation is less clear for the chiral transition, although also here we observe a steeper fall for 
the SU{3)-case. The behavior of the quark condensate for temperatures below the critical one, i.e. 
the rise with temperature combined with the sharp drop at has also been seen in quenched lattice 
calculations [23]. Nevertheless, it may very well be that the quantitative aspects of this rise are subject 
to the systematic uncertainties of the lattice gluon data [20]. These uncertainties are also reflected in 
the 'noisy' behavior of the quark condensate. Nevertheless it is remarkable that below Tc the dressed 
Polyakov loop is consistent with zero, signaling conserved centre symmetry. 

4.2 Unquenched QCD at finite temperature and chemical potential 

We now include two flavors of quarks via the quark loop as explained above. The effect of the matter 
sector on the gluon is a reduction of the dressing functions, which leads to a reduced interaction strength 
in the quark self-energy, and therefore to a smaller critical temperature. In Fig. H] we show the evolution 
of the order parameters at /i = 0. What we find is a crossover for both the condensate and the dressed 
Polyakov loop. The value for the pseudo-critical temperature is Tc ^~ = 180 ± 5 MeV from the quark 
condensate and Tc ^~ = 195 ±5 MeV from the dressed Polyakov loop. The difference in these numbers 
can be attributed to the crossover nature of the transition. 

When we go to /i > the condensate for neither periodic {ip = 0) nor anti-periodic (v? = tt) boundary 
conditions becomes complex. This leads to a difference in S+i and i.e. in the dressed Polyakov 
loop and its conjugate. In Fig. |5]the resulting phase diagram of two flavor QCD is shown. For the chiral 
transition we observe a crossover up to relatively large values of the chemical potential where we find a 
critical endpoint at (TEp,fiEp) ~ (95,280) MeV. Since i^ep/Tep ~ 3 1, this result suggests that the 
CEP is outside the reach of lattice QCD. For the confinement /deconfinement transition we observe that 
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Figure 5: The phase boundary for chiral symmetry and confinement at real chemical potential. The 
solid lines above the CEP denote the spinodals which mark the area of coexistence of chiral symmetric 
and broken solutions of the DSE. 

the critical temperature extracted from the dressed Polyakov loop and its conjugate is nearly equal, and 
close to that extracted from the quark condensate. As the chemical potential is increased the crossover 
becomes steeper and the two transition lines come closer together, meeting at around /i ~ 200 MeV. 

Both results, the CEP at large /x and the coinciding phase transitions agree well with results from 
the PQM model [3] beyond mean field, where the matter back-reaction on the Yang-Mills sector is also 
taken into account. 

We should note here that at the chemical potentials where we find the CEP our truncation scheme 
becomes less reliable, since the infiuence of Baryons is neglected. It may therefore be advised to rephrase 
our results as an exclusion of the CEP in the /i/T < 1 region. This is consistent with longstanding 
predictions from investigations of the curvature of the chiral critical surface in the Columbia plot [21] 
and also with recent lattice results on the curvature of the chiral and deconfinement crossover lines at 
small chemical potential [25]. 

5 Conclusion 

We have presented a truncation scheme for the Dyson- Schwinger equations of QCD where we take data 
from a lattice calculation for the temperature dependent quenched gluon, and introduce the quark loop 
for studies of unquenched QCD, namely at finite chemical potential. 

Within this truncation we investigated the behavior of the quark condensate as an order parameter 
for chiral symmetry breaking, and of the dressed Polyakov loop as an order parameter for confinement. 
In the quenched case at /i = we found that the order parameters reproduce the lattice input, hinting 
at a second order phase transition for SU{2) and a first order phase transition for SU{3). At finite 
density we found that thermal fiuctuations from the matter sector lead to a critical endpoint at large 
densities while chiral and deconfinement transitions coincide. Our results serve as a basis for further 
studies of hot and dense QCD. 
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